Genome-wide host-pathogen analyses reveal genetic interaction points in tuberculosis disease.

The genetics underlying tuberculosis (TB) pathophysiology are poorly understood. Human genome-wide association studies have failed so far to reveal reproducible susceptibility loci, attributed in part to the influence of the underlying Mycobacterium tuberculosis (Mtb) bacterial genotype on the outcome of the infection. Several studies have found associations of human genetic polymorphisms with Mtb phylo-lineages, but studies analysing genome-genome interactions are needed. By implementing a phylogenetic tree-based Mtb-to-human analysis for 714 TB patients from Thailand, we identify eight putative genetic interaction points (P < 5 × 10-8) including human loci DAP and RIMS3, both linked to the IFNγ cytokine and host immune system, as well as FSTL5, previously associated with susceptibility to TB. Many of the corresponding Mtb markers are lineage specific. The genome-to-genome analysis reveals a complex interactome picture, supports host-pathogen adaptation and co-evolution in TB, and has potential applications to large-scale studies across many TB endemic populations matched for host-pathogen genomic diversity.

Transferability of genetic loci and polygenic scores for cardiometabolic traits in British Pakistani and Bangladeshi individuals.

Individuals with South Asian ancestry have a higher risk of heart disease than other groups but have been largely excluded from genetic research. Using data from 22,000 British Pakistani and Bangladeshi individuals with linked electronic health records from the Genes & Health cohort, we conducted genome-wide association studies of coronary artery disease and its key risk factors. Using power-adjusted transferability ratios, we found evidence for transferability for the majority of cardiometabolic loci powered to replicate. The performance of polygenic scores was high for lipids and blood pressure, but lower for BMI and coronary artery disease. Adding a polygenic score for coronary artery disease to clinical risk factors showed significant improvement in reclassification. In Mendelian randomisation using transferable loci as instruments, our findings were consistent with results in European-ancestry individuals. Taken together, trait-specific transferability of trait loci between populations is an important consideration with implications for risk prediction and causal inference.

Low Levels of Factor H Family Proteins During Meningococcal Disease Indicate Systemic Processes Rather Than Specific Depletion by Neisseria meningitidis.

Neisseria meningitidis, the causative agent of meningococcal disease (MD), evades complement-mediated clearance upon infection by 'hijacking' the human complement regulator factor H (FH). The FH protein family also comprises the homologous FH-related (FHR) proteins, hypothesized to act as antagonists of FH, and FHR-3 has recently been implicated to play a major role in MD susceptibility. Here, we show that the circulating levels of all FH family proteins, not only FH and FHR-3, are equally decreased during the acute illness. We did neither observe specific consumption of FH or FHR-3 by N. meningitidis, nor of any of the other FH family proteins, suggesting that the globally reduced levels are due to systemic processes including dilution by fluid administration upon admission and vascular leakage. MD severity associated predominantly with a loss of FH rather than FHRs. Additionally, low FH levels associated with renal failure, suggesting insufficient protection of host tissue by the active protection by the FH protein family, which is reminiscent of reduced FH activity in hemolytic uremic syndrome. Retaining higher levels of FH may thus limit tissue injury during MD.

Systematic analysis of Kaposi's sarcoma (KS)-associated herpesvirus genomes from a KS case-control study in Cameroon: Evidence of dual infections but no association between viral sequence variation and KS risk.

In sub-Saharan Africa, Kaposi's sarcoma-associated herpesvirus (KSHV) is endemic, and Kaposi's sarcoma (KS) is a significant public health problem. Until recently, KSHV genotype analysis was performed using variable gene regions, representing a small fraction of the genome, and thus the contribution of sequence variation to viral transmission or pathogenesis are understudied. We performed near full-length KSHV genome sequence analysis on samples from 43 individuals selected from a large Cameroonian KS case-control study. KSHV genomes were obtained from 21 KS patients and 22 control participants. Phylogenetic analysis of the K1 region indicated the majority of sequences were A5 or B1 subtypes and all three K15 alleles were represented. Unique polymorphisms in the KSHV genome were observed including large gene deletions. We found evidence of multiple distinct KSHV genotypes in three individuals. Additionally, our analyses indicate that recombination is prevalent suggesting that multiple KSHV infections may not be uncommon overall. Most importantly, a detailed analysis of KSHV genomes from KS patients and control participants did not find a correlation between viral sequence variations and disease. Our study is the first to systematically compare near full-length KSHV genome sequences between KS cases and controls in the same endemic region to identify possible sequence variations associated with disease risk.

Integrating polygenic risk scores in the prediction of type 2 diabetes risk and subtypes in British Pakistanis and Bangladeshis: A population-based cohort study.

BACKGROUND: Type 2 diabetes (T2D) is highly prevalent in British South Asians, yet they are underrepresented in research. Genes & Health (G&H) is a large, population study of British Pakistanis and Bangladeshis (BPB) comprising genomic and routine health data. We assessed the extent to which genetic risk for T2D is shared between BPB and European populations (EUR). We then investigated whether the integration of a polygenic risk score (PRS) for T2D with an existing risk tool (QDiabetes) could improve prediction of incident disease and the characterisation of disease subtypes. METHODS AND FINDINGS: In this observational cohort study, we assessed whether common genetic loci associated with T2D in EUR individuals were replicated in 22,490 BPB individuals in G&H. We replicated fewer loci in G&H (n = 76/338, 22%) than would be expected given power if all EUR-ascertained loci were transferable (n = 101, 30%; p = 0.001). Of the 27 transferable loci that were powered to interrogate this, only 9 showed evidence of shared causal variants. We constructed a T2D PRS and combined it with a clinical risk instrument (QDiabetes) in a novel, integrated risk tool (IRT) to assess risk of incident diabetes. To assess model performance, we compared categorical net reclassification index (NRI) versus QDiabetes alone. In 13,648 patients free from T2D followed up for 10 years, NRI was 3.2% for IRT versus QDiabetes (95% confidence interval (CI): 2.0% to 4.4%). IRT performed best in reclassification of individuals aged less than 40 years deemed low risk by QDiabetes alone (NRI 5.6%, 95% CI 3.6% to 7.6%), who tended to be free from comorbidities and slim. After adjustment for QDiabetes score, PRS was independently associated with progression to T2D after gestational diabetes (hazard ratio (HR) per SD of PRS 1.23, 95% CI 1.05 to 1.42, p = 0.028). Using cluster analysis of clinical features at diabetes diagnosis, we replicated previously reported disease subgroups, including Mild Age-Related, Mild Obesity-related, and Insulin-Resistant Diabetes, and showed that PRS distribution differs between subgroups (p = 0.002). Integrating PRS in this cluster analysis revealed a Probable Severe Insulin Deficient Diabetes (pSIDD) subgroup, despite the absence of clinical measures of insulin secretion or resistance. We also observed differences in rates of progression to micro- and macrovascular complications between subgroups after adjustment for confounders. Study limitations include the absence of an external replication cohort and the potential biases arising from missing or incorrect routine health data. CONCLUSIONS: Our analysis of the transferability of T2D loci between EUR and BPB indicates the need for larger, multiancestry studies to better characterise the genetic contribution to disease and its varied aetiology. We show that a T2D PRS optimised for this high-risk BPB population has potential clinical application in BPB, improving the identification of T2D risk (especially in the young) on top of an established clinical risk algorithm and aiding identification of subgroups at diagnosis, which may help future efforts to stratify care and treatment of the disease.

The genomics of heart failure: design and rationale of the HERMES consortium.

AIMS: The HERMES (HEart failure Molecular Epidemiology for Therapeutic targetS) consortium aims to identify the genomic and molecular basis of heart failure. METHODS AND RESULTS: The consortium currently includes 51 studies from 11 countries, including 68 157 heart failure cases and 949 888 controls, with data on heart failure events and prognosis. All studies collected biological samples and performed genome-wide genotyping of common genetic variants. The enrolment of subjects into participating studies ranged from 1948 to the present day, and the median follow-up following heart failure diagnosis ranged from 2 to 116 months. Forty-nine of 51 individual studies enrolled participants of both sexes; in these studies, participants with heart failure were predominantly male (34-90%). The mean age at diagnosis or ascertainment across all studies ranged from 54 to 84 years. Based on the aggregate sample, we estimated 80% power to genetic variant associations with risk of heart failure with an odds ratio of ≥1.10 for common variants (allele frequency ≥ 0.05) and ≥1.20 for low-frequency variants (allele frequency 0.01-0.05) at P < 5 × 10-8 under an additive genetic model. CONCLUSIONS: HERMES is a global collaboration aiming to (i) identify the genetic determinants of heart failure; (ii) generate insights into the causal pathways leading to heart failure and enable genetic approaches to target prioritization; and (iii) develop genomic tools for disease stratification and risk prediction.

Identification of novel locus associated with coronary artery aneurysms and validation of loci for susceptibility to Kawasaki disease.

Kawasaki disease (KD) is a paediatric vasculitis associated with coronary artery aneurysms (CAA). Genetic variants influencing susceptibility to KD have been previously identified, but no risk alleles have been validated that influence CAA formation. We conducted a genome-wide association study (GWAS) for CAA in KD patients of European descent with 200 cases and 276 controls. A second GWAS for susceptibility pooled KD cases with healthy paediatric controls from vaccine trials in the UK (n = 1609). Logistic regression mixed models were used for both GWASs. The susceptibility GWAS was meta-analysed with 400 KD cases and 6101 controls from a previous European GWAS, these results were further meta-analysed with Japanese GWASs at two putative loci. The CAA GWAS identified an intergenic region of chromosome 20q13 with multiple SNVs showing genome-wide significance. The risk allele of the most associated SNV (rs6017006) was present in 13% of cases and 4% of controls; in East Asian 1000 Genomes data, the allele was absent or rare. Susceptibility GWAS with meta-analysis with previously published European data identified two previously associated loci (ITPKC and FCGR2A). Further meta-analysis with Japanese GWAS summary data from the CASP3 and FAM167A genomic regions validated these loci in Europeans showing consistent effects of the top SNVs in both populations. We identified a novel locus for CAA in KD patients of European descent. The results suggest that different genes determine susceptibility to KD and development of CAA and future work should focus on the function of the intergenic region on chromosome 20q13.

Distinct genetic architectures and environmental factors associate with host response to the γ2-herpesvirus infections.

Kaposi's sarcoma-associated herpesvirus (KSHV) and Epstein-Barr Virus (EBV) establish life-long infections and are associated with malignancies. Striking geographic variation in incidence and the fact that virus alone is insufficient to cause disease, suggests other co-factors are involved. Here we present epidemiological analysis and genome-wide association study (GWAS) in 4365 individuals from an African population cohort, to assess the influence of host genetic and non-genetic factors on virus antibody responses. EBV/KSHV co-infection (OR = 5.71(1.58-7.12)), HIV positivity (OR = 2.22(1.32-3.73)) and living in a more rural area (OR = 1.38(1.01-1.89)) are strongly associated with immunogenicity. GWAS reveals associations with KSHV antibody response in the HLA-B/C region (p = 6.64 × 10-09). For EBV, associations are identified for VCA (rs71542439, p = 1.15 × 10-12). Human leucocyte antigen (HLA) and trans-ancestry fine-mapping substantiate that distinct variants in HLA-DQA1 (p = 5.24 × 10-44) are driving associations for EBNA-1 in Africa. This study highlights complex interactions between KSHV and EBV, in addition to distinct genetic architectures resulting in important differences in pathogenesis and transmission.

Genome-Wide Sequence Analysis of Kaposi Sarcoma-Associated Herpesvirus Shows Diversification Driven by Recombination.

Background: Kaposi sarcoma-associated herpesvirus (KSHV) establishes lifelong infection in the human host and has been associated with a variety of malignancies. KSHV displays striking geographic variation in prevalence, which is highest in sub-Saharan Africa. The current KSHV genome sequences available are all tumor cell line-derived or primary tumor-associated viruses, which have provided valuable insights into KSHV genetic diversity. Methods: Here, we sequenced 45 KSHV genomes from a Ugandan population cohort in which KSHV is endemic; these are the only genome sequences obtained from nondiseased individuals and of KSHV DNA isolated from saliva. Results: Population structure analysis, along with the 25 published genome sequences from other parts of the world, showed whole-genome variation, separating sequences and variation within the central genome contributing to clustering of genomes by geography. We reveal new evidence for the presence of intragenic recombination and multiple recombination events contributing to the divergence of genomes into at least 5 distinct types. Discussion: This study shows that large-scale genome-wide sequencing from clinical and epidemiological samples is necessary to capture the full extent of genetic diversity of KSHV, including recombination, and provides evidence to suggest a revision of KSHV genotype nomenclature.

High genotypic diversity among rotavirus strains infecting Gambian children.

BACKGROUND: Rotavirus is the leading cause of diarrhea in children <5 years of age. In light of the implementation of rotavirus vaccines of limited valency, it is important to characterize the genotypic diversity of circulating rotavirus in sub-Saharan Africa. METHODS: We collected stool samples from children 0-59 months of age who presented at the health centres as cases with moderate-to-severe diarrhea in the Upper River Region of The Gambia. Stool samples were also collected from age, sex and area-matched healthy controls. All stool samples were assayed for rotavirus antigens by enzyme-linked immunosorbent assay and genotyping was done using reverse transcriptase polymerase chain reaction. RESULTS: We enrolled 1029 cases and 1569 controls during the 3-year study period (2008-2010). The detection rate of rotavirus among the cases was 20% (204/1029) and 3% (42/1569) among controls. At least 18 genotypes were found and the predominant genotypes were G2P[6] (28%), G1P[8] (26%) and G1P[10] (10%). The rare identified genotypes (<1%) were G2P[14], G8P[6], G9P[6] and G4P[10]. There was also a strong positive association between rotavirus infection and the dry season (odds ratio: 9.83, 95% confidence interval: 6.18-15.63, P < 0.001). A significant increase in the odds of rotavirus and G1P[8] detection with the use of untreated water and the presence of cats, rodents and cows in the child's residence was also found. CONCLUSION: This study provides important baseline data for the genotypes circulating before vaccine implementation. The wide diversity of genotypes circulating in The Gambia implies the need for vigilant effectiveness surveillance following the implementation of RotaTeq in August 2013.

Immunogenic Mycobacterium africanum strains associated with ongoing transmission in The Gambia.

In West Africa, Mycobacterium tuberculosis strains co-circulate with M. africanum, and both pathogens cause pulmonary tuberculosis in humans. Given recent findings that M. tuberculosis T-cell epitopes are hyperconserved, we hypothesized that more immunogenic strains have increased capacity to spread within the human host population. We investigated the relationship between the composition of the mycobacterial population in The Gambia, as measured by spoligotype analysis, and the immunogenicity of these strains as measured by purified protein derivative-induced interferon-γ release in ELISPOT assays of peripheral blood mononuclear cells. We found a positive correlation between strains with superior spreading capacity and their relative immunogenicity. Although our observation is true for M. tuberculosis and M. africanum strains, the association was especially pronounced in 1 M. africanum sublineage, characterized by spoligotype shared international type 181, which is responsible for 20% of all tuberculosis cases in the region and therefore poses a major public health threat in The Gambia.

Chrysomya putoria, a putative vector of diarrheal diseases.

BACKGROUND: Chrysomya spp are common blowflies in Africa, Asia and parts of South America and some species can reproduce in prodigious numbers in pit latrines. Because of their strong association with human feces and their synanthropic nature, we examined whether these flies are likely to be vectors of diarrheal pathogens. METHODOLOGY/PRINCIPAL FINDINGS: Flies were sampled using exit traps placed over the drop holes of latrines in Gambian villages. Odor-baited fly traps were used to determine the relative attractiveness of different breeding and feeding media. The presence of bacteria on flies was confirmed by culture and bacterial DNA identified using PCR. A median of 7.00 flies/latrine/day (IQR = 0.0-25.25) was collected, of which 95% were Chrysomya spp, and of these nearly all were Chrysomya putoria (99%). More flies were collected from traps with feces from young children (median = 3.0, IQR = 1.75-10.75) and dogs (median = 1.50, IQR = 0.0-13.25) than from herbivores (median = 0.0, IQR = 0.0-0.0; goat, horse, cow and calf; p<0.001). Flies were strongly attracted to raw meat (median = 44.5, IQR = 26.25-143.00) compared with fish (median = 0.0, IQR = 0.0-19.75, ns), cooked and uncooked rice, and mangoes (median = 0.0, IQR = 0.0-0.0; p<0.001). Escherichia coli were cultured from the surface of 21% (15/72 agar plates) of Chrysomya spp and 10% of these were enterotoxigenic. Enteroaggregative E. coli were identified by PCR in 2% of homogenized Chrysomya spp, Shigella spp in 1.4% and Salmonella spp in 0.6% of samples. CONCLUSIONS/SIGNIFICANCE: The large numbers of C. putoria that can emerge from pit latrines, the presence of enteric pathogens on flies, and their strong attraction to raw meat and fish suggests these flies may be common vectors of diarrheal diseases in Africa.